A van der Corput-type algorithm for LS-sequences of points 

Ingrid CarboneQ 

Abstract 

In this paper we associate to any LS-sequence of partitions {p'l s } the 
corresponding LS-sequence of points {£" s } obtained reordering the points 
of each partition with an explicit algorithm. The procedure begins with the 
representation in base L + S of natural numbers, [h]l+s. and ends with the LS- 
radical inverse function <j>i t s, introduced ad hoc, evaluated at an appropriate 
subsequence of natural numbers depending on L and 5. This construction 
is deeply related to the geometric representation of the points of by 
suitable affine functions and reminds the van der Corput sequences in base 
b. 
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1 Introduction 

The search for new classes of uniformly distributed sequences of points on [0, 1 [ or 
in higher dimension is one of the most interesting aspects of the theory of uniform 
distribution, especially after the discovery that low discrepancy sequences can be 
applied to Quasi Monte-Carlo Methods. Some of them are naturally associated to 
sequences of partitions. The most important one is the van der Corput sequence, 
introduced in lfT6ll . which is a low discrepancy sequence of points obtained by 
choosing in the "right" order the left endpoints of the binary sequence of partitions 
{8 n }, where for each n G N 

,0<i<2" 1 j. 

We will come back to van der Corput sequence very soon. 

In a recent paper [2] the author has introduced a new countable family of uni- 
formly distributed sequences of partitions, called LS-sequences of partitions and 
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denoted by {p1 s }\ in the simple case L = S = 1 we get one of the Kakutani se- 
quences of partitions; we have called it Kakutani-Fibonacci sequence of partitions 
and we will describe it in a very deatailed way in Section 3 and Section 4. Upper 
and lower estimates of the discrepancy of {p'l s } have been also given and it has 
been showed that {p'l s } has low discrepancy whenever L > S; in particular, the 
Kakutani-Fibonacci sequence has low discrepancy (it is, actually, the first time 
that a precise estimate of the discrepancy of one of the Kakutani sequences has 
been obtained). 

By reordering the left endpoints of the intervals of each partition p£ s with an 
explicit algorithm described in the same article, it is has been possible to asso- 
ciate to each sequence of partition {p^ s } a sequence of points called LS-sequence 
of points and denoted by Upper bounds for the discrepancy of these se- 

quences have been given, and the proof that each {^5} is a low discrepancy 
sequence when L> S, too. More details on {pl s }, {^5} an d their discrepancy 
will be given in Section 2. 

In Section 3 we present a new algorithm to construct LS- sequences of points 
{^l s }, which is formally independent at all on the concept of sequence of par- 
titions: we do not need to reorder the points of each partitions, as it is shown in 
Theorem 13.51 and therefore it is simpler to use then the one introduced in 0. 
Of course, the connection between each {^[ s } and the corresponding {p'l s } lies 
behind all the procedure, and this geometric interpretation represents one of the 
most interesting aspects of this new family of sequences. 

In order to prove this result, we represent the points of each partition p£ s as 
subsequent compositions of suitable functions depending on L and S, as showed 
in Proposition [X2j 

We conclude this paper with a significative example. 

Let us complete this Introduction with some classical definitions, starting from 
the constructions of the van der Corput sequence. 
Any natural number n can be expressed as 

M 

"=£4)2*, (1) 

k=0 

with a k (n) E {0, 1} for all < k < M, with M = [log 2 n] (here and in the sequel [•] 
denotes the integer part). Such representation is unique modulus the addition of 
higher powers of 2 with zero coefficients. 

The expression (OQ) leads to the binary representation 

[n}2 = a M (n)aM-i(n)...ao(n). (2) 
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The representation of n in base 2 given by (OQ) is used to define the radical- 
inverse function 02 on N which associates to the string of digits © the number 

M 

<h(n) = y £a k (n)2- k -\ 

k=Q 

whose binary representation is 0.aQ{n)a\(n) ...aM(n). 

Of course < tyiin) < 1 for all n > and the 2-radix notation of 02 (ft) actually 
is [(fe(n)]2 = 0.ao(n)ai(n)...aM(n). 

The van der Corput sequence is the sequence 



{^Wjn- 
After its introduction in 1935, several generalizations of the van der Corput se- 
quence have been given during the subsequent decades. Among them, we mention 
the so called van der Corput sequence in base b denoted by {(j)b(n)} n , where b is 
any positive integer, with <fc(n) = Y!k=Q^k{n)b ~ k _1 and ai[n) E {0, l,...,b— 1}. 
These sequences were introduced by Hammersley ( ifTTTO and their discrepancy 
was first studied by Halton ( ifTOlO . who proved that these sequences have low dis- 
crepancy for all b. 

Another generalization is the so called scrambled van der Corput sequence 
denoted by {§°{n)}„, where $°{n) = Lf=o G k( a k( n ))b~ k ~ l and {a n } n is a se- 
quence of permutations of the set {0, 1, . . . ,Z? — 1}. Faure ((El and flU) proved 
that, for some suitable {(7 n }„, the discrepancy of {§°{n)} n is asymptotically bet- 
ter that the discrepancy of {§b(n)} n . 

We need some definitions (see and lfT4lO . 

Definition 1.1. Given a sequence of points {x n } of the interval [0, 1[, we say that 

{x n } is uniformly distributed ( u.d.) if for any a, b E R such that < a < b < lwe 
have 

1 N 

Definition 1.2. Given a point set Xn = (xi,X2, ... ,xjv) of the interval [0, 1 [, the 

discrepancy of Xn is 



D(X N ) = sup 

Q<a<b<\ 



1 N 

^JLX[a,b[(Xj)-(b- 
ly 7=1 
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and the star-discrepancy of is 



D*(X N ) 



sup 

0<b<l 



1 N 

xLX[o,b[( x j)- b 



Given a sequence of points X = {x n } in [0, 1[, for each N we can consider 
the point set = (x\,X2, ■ ■ ■ ,xjv) and the associated sequences {D(Xn)}n (and, 
respectively, {D*(Xjv)}/v)- It is important to study how these sequences behave 
when N tends to infinity. 

It is well known that a sequence X = {x„} is u.d. if and only if the discrep- 
ancy D{Xn) tends to whenever N tends to infinity, and also that there exists an 
absolute constant c such that c logAf < ND{X N ) < N (see lfT5l0 . 

We say that X = {x n } has low discrepancy if there exists a constant C > such 
that ND(X N ) < C log N for any N. 

If we consider a sequence of partitions {K n } of [0, 1 [, with 

% = {[yf,yW 1 [,l<7<a ) 

each partition % n can be regarded as a point set. For this reason, we can re- write the 
definition of uniform distribution and discrepancy, given for sequences of points, 
in this particular setting. 

We say that the sequence of partitions is u.d. if for any real numbers a, b 
such that 0<a<b<lwe have 



l n j—i 



a. 



Moreover, since each % n is defined by t n points, the concepts of discrepancy 
and star-discrepancy can be naturally extended to the sequence {n n } as follows: 



D(n n 



sup 

0<a<b<l 



rLxia, b[ (yj )-(*>■ 

ln 7=1 



and 



sup 

0<b<l 



We say that {n n } has low discrepancy if there exists a constant C > such that 
t n D(n n ) < C for any n (see 0). 



4 



In IfTTl it has been introduced a new class of u.d. sequences of partitions called 
p -refinements. 

Definition 1.3. For any fixed non trivial finite partition p of [0, 1 [, the p -refinement 
of a partition % of [0, 1 [ (denoted by p%) is obtained by subdividing all the intervals 
of % having maximal length positively (or directly) homothetically to p. If for 
any «eNwe denote by p n Tt the p-refinement of p n ~ l n, we get a sequences of 
partitions {p n 7t}, called the sequence of successive p-refinements of n. 

If p = {[0, a[, [a, 1[} for some a e [0, 1[ and CO = {[0, 1[} is the trivial par- 
tition, the sequence {p n O)} is actually the well known Kakutani's sequence of 
a-refinements, which has been defined in ifPIll and will be denoted by { K n } . Kaku- 
tani proved that 

Theorem 1.4. The sequence {K n } is uniformly distributed. 

In IfTTl it has been proved the following generalization: 
Theorem 1.5. The sequence {p n (o} is uniformly distributed. 

The uniform distribution of the sequence {p n n}, when % is not the trivial par- 
tition, has been investigated in 0]|, where nice necessary and sufficient conditions 
on % and p have been found. 

A further contribution to the theory of p-refinement, moreover, has been given 
by Drmota and Infusino who, in J6l, provided estimates of the discrepancy of 
sequences of p-refinements and, in particular, of LS- sequences of partions. 

In the last years, however, the interest on Kakutani's splitting procedure gave 
birth to several generalizations, among which we would like to mention the ex- 
tension to separable metric spaces [3] and the generalization to higher dimensions 
which is intrinsically ^-dimensional [4J. Moreover, it has been extended to some 
fractals in lfl2l 

In [[5]| the authors present a von Neumann-type theorem for sequences of parti- 
tions, showing actually that a sequence of partitions having an infinitesimal diam- 
eter can be reordered in a sequence of partitions which is uniformly distributed. 

2 LS-sequences 

In this section we recall the definition of LS-sequence of points {££ s } introduced 
in 0. To this aim, we need the definition of LS-sequence of partitions {p'[ s }- 
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Definition 2.1. Let us fix two positive integers L and S and let 7 be the positive 
solution of Ly+Sy 2 = 1 . Denote by pi 5 the partition defined by L "long" intervals 
having length 7 followed by S "short" intervals having length y 2 . The sequence 
of the successive ^-refinements of the trivial partition CO is called LS-sequence 
of partitions and is denoted by {p'l s co} (or {pl s } for short). 

If we denote with t„ the total number of intervals of p£ s , with /„ the number 
of its long intervals and with s n the number of its short intervals, it is very simple 
to see that t n = l„ + s n , l n = Ll n -\ + s n -i and s n = Sl n -\. We deduce that t n = 
Lt n -i +St n -2 (with to = 1 and t\ = L + S) which has the explicit solution 



l+5 7 /IV Sy-Sf 



We note that l n and s n satisfy the same difference equation, with initial conditions 
Iq = 1, l\ = L and so = 0,si = S, respectively. 

In |[2l we have given precise estimates from above and from below of the 
discrepancy of LS-sequences of partitions, and proved the following 

Theorem 2.2. i)IfS<L there exist ci,c[ > such that for any n e N 

c[<t n D(pl s )< Cl . 



ii) IfS = L+l there exist C2, c' 2 > such that for any neN 

c' 2 logt n < t n D(pl s ) < c 2 \op n . 



Hi) IfS>L + 2 there exist C3, C3 > such that for any «eN 

cW" <tnD(pl s ) < c 3 t l n -\ 

where l-r = - 1 ^f>0. 

We emphasize that each L5-sequence of partitions with S < L has low discrep- 
ancy. 

Let us now recall the definition of the sequence of points {^5} associated to 
the sequence of partitions {p'l s }- 
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Definition 2.3. We fix a sequence of partitions {p'l s }. For each n>lwe define 
the two families of functions 

( P l n+1 \x)=x + if+ 1 and <$j l \x) = x + Lf +1 +jf+ 2 (3) 

for 1 < i < L and 1 < j < S — 1. We denote by A| s the set made by the t\ left 

endpoints of the intervals of p[ s ordered by magnitude. If A n L s = (j;^ , . . . , <^f^ 

denotes the set of the t n points defining pi s , ordered in the appropriate way, the 

set A£~t , which consists in a reordering of the points defining p£ ^ , is obtained 
as follows: 

\n+l _ (An) An) An) 

ft ^ i »(tf»),...,fr i \« k w ),..,ft^ ,) (« 1 w ),..,ft M (fif > ). 
f^ l, «. w ).- • • • fir'Kf). .... #.\ > «. w >- • • , &t>)) • (4> 

The sequence obtained reordering successively the points of p'l + s l \p£s * s called 
LS-sequence of points and is denoted by {^£5}- 

In we also provided the following bounds for the discrepancy of these 
sequences of points. 

Theorem 2.4. i) IfS < L there exists k\>§ such that for any WgN 

ii) IfS = L+l there exists k2, c' 2 > such that for any iV6N 

c' 2 \ogN<ND (<^ s , ^ ■ ■ ■ , ft) < ^2 log 2 AT . 
Hi) IfS>L + 2 there exists £3, c 3 > smc/i that for any N G N 

w/i e r e l-T = -i^>0. 
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Theorem 12.21 and Theorem 12.41 show that if the LS-sequence of partitions has 
low discrepancy, the corresponding LS-sequence of points has low discrepancy, 
too. Moreover, c' 2 and c' 3 denote the same constants in both theorem as, of course, 

whenever N = t„ for some n, it follows that D s , s , . . . , %ls) = ^(Pl s) ■ 

Note that the lower bounds just follow from Theorem 12.21 In fact we believe 
that the discrepancy of all the LS- sequences of points coincides with the upper 
bound given by Theorem [231 

3 The new definition " a la van der Corput" 

From now on we fix L, S and 7 such that Ly+Sy 2 = 1 . 

The first step is to represent the points of the sequence {^5} writing the 
elements of each set A n L s as images of the point x = under the compositions of 
suitable functions we are going to define. 

To this purpose we introduce the functions y/i on restricted domains as follows: 
for every < i < L — 1 we consider 

Wi{x) = Jx + iy restricted to < x< 1, (5) 

while for every L<i<L+S— lwe define 

Wi (x) = yx J rLy J r{i-L)y 2 restricted to < x < y. (6) 

We observe that the functions © map [0, 1 [ onto [iy, (i + l)y [ and the func- 
tions © map [0, y[onto [Ly+(i— L)y 2 1 Ly+(i— L+ l)y 2 [. Therefore, the compo- 
sitions \jfi o y/j are not defined if and only if L<i < L+S — 1 and 1 < j < L+S— 1. 

Let us denote by Ei s the set consisting of all the pairs of indices which corre- 
spond to the "forbidden" compositions, i.e. 

£ L ,s = {L,L+l,...,L + S-l}x{l,...,L + S-l}. (7) 

In order to present our main result, we have to give the explicit expression of 
the compositions yi l i 2 ...i„ = Wh tyi 2 ° ■■■ Win °f me functions © and ©. 

Lemma 3.1. For any n G N and any n-tuple (hi2---in) of elements of the set 
{0, 1, . . . ,L+S— 1} such that (ihJh+i) ^ for any 1 <h<n—\, we have 

Vnn-M = fx+thr k , (8) 
k=\ 

where = it if < < L— 1, while = L+ — L)y if L < i\ < L+ S — 1. 
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Proof. We prove © by induction on n. 

If n = 2, because of © and © we have only three kinds of compositions 

a) if < i\ , 12 < L — 1, we have ijf^ j 2 (x) = y 2 x + £^ =1 y*, where b^ = i^; 

b) if < i\ < L — 1 and L<i2<L + S — 1, a simple calculation gives ®, with 
= i] and b2 = L+(i2 — L) y; 

c) if L < z',- < L + S — 1 and < z'2 < L — 1 , we obtain the same expression, where 
b\ = L + (z'i — L)y and b2 = 12- 

Suppose now that (fTTT) holds for n > 1 and let us prove it for n + 1, namely 

that 

n+l 

vw..%;„ + iM = f + ^I>y* (9) 

for all z, 1+ i 7^ 0. Of course, at this stage we have to distinguish between the two 
kinds of functions © and ©. 

If < i n+ i < L— 1, we simply obtain © with b n+ \ = i n+ {. 

If L < i n+ i < L + S — 1, simple calculations give again ©, with b n+ i = L + 
(i n +i —L)y. So the induction is complete. □ 

An immediate consequence of the previous lemma is the following result. 
Proposition 3.2. The first t n points of the LS-sequence {^l s } are 

* n L ,s=Whh...i n m, do) 

where (ihJh+i) ^ £l,s for any 1 < h < n — 1 and (i n , i n -i, - ■ ■ , ii) are the n-tuples 
of elements of the set {0, 1 , . . . , L + S — 1 }, ordered with respect to the magnitude 
of the numbers i n i n - 1 • • • h in base L + S. 

Proof. Taking © and © into account, we see that 

Ai, s = {Vb(0),Vi(0),.--,VL+s-i(0)}, 

which are exactly the L + S left endpoints of the intervals of p[ s ordered by mag- 
nitude, so (flOl ) is true for n = 1. 

Fix any integer n > 1 and suppose (flOl ) is true for n. 

If i n+ \ = 0, as i//o(0) = we have 
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and we simply observe that this way we get all the points of A n LS . All the new 

points of A n L + l \A n L s are obtained when i n+ \ ^ as follows. 
If 1 < i n+ \ < L we have 

n 

Yhi 2 ...i n i n+1 (0) = i n +if +1 + = VniiMty+in+if^ 

according to @. 

If L + 1 < i n+ i < L + S — 1, we set j n+ \ = i n+ \ —L, with 1 < j n+ \ < S— 1, and 
write 

Vm 2 ..,„wi(°) = (L+j n+ ir)f +1 + thy k 

k=i 

= xif h i,..i„(0) +Lf +l +j n+ ,f +2 = cp^( ¥ni2 ... in (0))- 

At this point we recall that the set A^y \ A2 s is made by (L + S — 1 )l n points 
as we evaluate the L + S — I functions defined by © at the first /„ points of A'[ s 
(see ©). 

It remains only to prove that the total number of the points of A L + S \A'l s we 

add in the way described above are exactly (L + S — l)l n . We denote by the 
total number of n-tuples containing two consecutive digits of the set Ei $ defined 

by ©. Of course, we have = (L + S) n - t n . 

We observe, in fact, that the points of A£ ^ \ A'[ s come from L + S —1 blocks 
of all the (n + l)-tuples of elements of the set {0, 1, . . . ,L + S — 1}, which are 
(L + S) n , ordered by magnitude of the numbers in base L + S. Moreover, the 
total number of not admissible (n + 1) -tuples which are contained in these blocks 
is given by the difference between (L + S— l)(L + S) n and the total number of 

forbidden ones, which is exactly af^' — d^/ s . In other words, 
(L + S-l)(L + S y-(d^ l) -d%) 

= (L + S-1)(L + S) n - ((L + S) n+1 -t n+ i -(L + S-l) n + t n ) 
= t n +\ —tn = (L + S— \)l n . 

The theorem is completely proved. □ 
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We are now ready to present LS-sequences of points from another point of 
view, which is deeply connected to the explicit geometric representations given 
in Proposition 13.21 More precisely, we want to show how to get the points of 
these sequences starting from the representation in base L + S of natural numbers. 
This algorithm is very simple and allows us to compute directly the points of each 

Given 

M 

n=Y,"k(n)(L+S) k , 

where a k (n) 6 {0, 1 , 2, . . . , L + S - 1 } for all < k < M and M = [log L+5 n] , its 
representation in base L + S is 

[n]L+s = ain(n) a M -i(n) ...a (n). 

If we reverse the order of the digits of [n}i + s an d a dd a dot, we get 

[h+s(n)] L+s = 0.a (n) a\ (n) ... a M -i(n) a M (n) , 

which is the representation in base L + S of the (L + S) -radical inverse function 
(j>L+s defined on N by 

M 

h+s(n) = Y,Mn)(L + S)- k - 1 . 

k=0 

Definition 3.3. We denote by Ni,s the set of all positive integers n, ordered by 
magnitude, with [h]l+s = «m(«)«m-i(«) •••«o(«) sucn that [a,k(n),ak+\{n)) £ 
E LjS foraU0<it<Af-l. 

Definition 3.4. For all n e N^s we define the LS-radical inverse function as fol- 
lows: 

M 

h,s( n ) = L d k(n)y k+ \ (11) 

k=0 

where dk(n) = akin) if < ajfc(n) < L—l and dk{n) = L + y(%(n) —L) if L < 
cik(n) < L + S — 1. 

Taking the above definitions into account, we are able to present the main 
result of this paper. 
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Theorem 3.5. Any LS-sequence of points coincides with the sequence {(j>L.s( n )} 
defined on Nl,s- 

Proof. Fix n e Nl,s- Taking Lemma |3~T1 and Definition 13 .41 into account, we write 

h,s( n ) = ¥n oni ...n M (0), 

which proves the theorem. □ 

Example 3.6. We consider the Kakutani-Fibonacci sequence of points l } cor- 
responding to L = S = 1 and y = j(\/5 — 1). See Fig. 1. 

1 6 4 3 8 2 7 5 

I 1 1 1 1 1 1 1 1 

y3 yl y2_|_y3y y_|_y4y_|_y3 \ 

Figure 1: The first 8 points of j}. 



o r i 



Figure 2: The functions ty/o and y/i associated to {^f x } 
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The functions defined by © and © reduce to y/o(x) = yx for < x < 1 and 
W\ (x) = yx + y for: < x < y. See the Fig. 2. 

According to Definition 13 .4[ Nij is the set of all natural numbers n such that 
the binary representation © does not contain two consecutive digits equal to 1. 
Moreover, the (1, 1) -radical inverse function defined by (PTTI) onN^i is 

M 

fo,i(n) = £ «Jt(*)/ +1 - 
k=o 

By Theorem [331 the Kakutani-Fibonacci sequence of points j} coincides 
with the sequence {0^1 (n)}. 

The following table shows the construction of the first 12 points of {^fj}, 
where the first column contains the first 12 elements of Nij. 



0^ 


[0]2 = 


->■ 


0.0 


[01,1 


(0)]2 


-> 





= a,i 


1 


[1]2 = 


1 -». 


0.1 


[01,1 


(1)]2 


->• 


7 


= ^i 


2^ 


[2]2 = 


10 -». 


0.01 


[01,1 


(2)]2 


->• 


r 




4^ 


[4]2 = 


100 -> 


0.001 = 


[01,1 


(4)] 2 


->• 


r 3 

7+7 3 




5^ 


[5]2 = 


101 ->■ 


0.101 = 


[01,1 


(5)]2 


— >■ 




8^ 


[8]2 = 


1000 ->■ 


0.0001 = 


[01,1 


(8)] 2 


->• 


r 4 




9^ 


[9] 2 = 


1001 ->■ 


0.1001 = 


[01,1 


(9)] 2 


— >■ 


7+/ 


= C 


10 


[10]2 = 


1010 -»■ 


0.0101 = 


[01,1 


(10)] 2 


->• 


y 2 + y 4 


= a 


16-> 


[16]2 = 


10000 -> 


0.00001 = 


[0i.i 


(16)] 2 


— >■ 


7 s 


- <E 9 
~~ 'i,i 


17-» 


m = 


10001 -> 


0.10001 = 


[01,1 


(17)]2 


->• 


y+y 5 


ElO 


18-> 


[18] 2 = 


10010 -». 


0.01001 = 


[0i,i 


(18)] 2 


-»• 


y 2 + y s 


- f 11 


20^ 


[20] 2 = 


10100 -> 


0.00101 = 


[01,1 


(20)] 2 


->• 


y 3 + y 5 


El2 

- Sl,l 
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